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We investigate light hadron spectroscopy with an improved quenched staggered quark 
action. We compare the results obtained with an improved gauge plus an improved 
quark action, an improved gauge plus standard quark action, and the standard gauge 
plus standard quark action. Most of the improvement in the spectroscopy results is due 
to the improved gauge sector. However, the improved quark action substantially reduces 
violations of Lorentz invariance, as evidenced by the meson dispersion relations. 



1 Introduction 

The precision of numerical lattice QCD simulations with the standard lattice actions is 
constrained by the available computational resources. In order to keep the duration of 
the calculation within manageable bounds, one is forced to use lattice spacings a which 
may be too large to accurately describe the continuum physics. This problem has been 
addressed by the development of improved |jl], |2|, Q and fixed-point 0, ^ actions. The 
promise of these actions is to yield good approximations to the continuum physics with 
relatively coarse lattice spacings. (For a recent review, see [|].) 

In the Symanzik improvement scheme |I|, Q the lattice action and fields are improved in 
powers of the lattice spacing a. This is achieved by introducing higher dimensional terms 
into the action. In the continuum limit, these terms are irrelevant, but at a finite cutoff 
the coefficients of these terms can be tuned so that the discretization errors of spectral 
quantities are diminished. The most straightforward method to determine the coefficients 
is to expand the action in a Taylor series in a, and cancel the leading scaling violating 
terms order by order (tree-level improvement). This can be refined by using perturbative 
analysis or non-perturbative numerical methods to determine the coefficients. 

In this paper, we study the improvement of the staggered (Kogut-Susskind) quark 
lattice QCD action with quenched spectroscopy calculations. The improvement is im- 
plemented by adding a third-nearest-neighbor term, first proposed by Naik more than a 
decade ago [0 . Some of the preliminary results of this study have already been published 
in 0, H . The same improvement scheme has been applied to nonzero temperature calcula- 



tions by Karsch et al. |T0[. The gauge configurations used in this study are generated with 
an O(a^) one-loop and tadpole-improved gauge action [^ |^. Since our main goal is to 
investigate the effects of the fermionic improvement, we compare the hadronic spectrum 
obtained with both the unimproved and improved fermion actions, using the same gauge 
configurations. An excellent baseline for the evaluation of the improvement is provided 
by our extensive standard (non-improved) quenched Kogut-Susskind hadron spectroscopy 
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As opposed to Wilson fermions, the improvement of the staggered action has attracted 
relatively little attention. This is partly due to the formal complexity of the staggered 
formulation, partly to the fact that the standard Wilson fermions have an error 0{a), 

1 



whereas the staggered action is already accurate to this order. Nevertheless, the improve- 
ment of the staggered action is highly desirable: the staggered action has a U(l)xU(l) 
chiral symmetry, remnant of the full continuum U(4)xU(4) symmetry (for 4 quark fla- 
vors). This symmetry is restored in the continuum limit; however, for practical values 
of the lattice spacing a substantial flavor symmetry breaking remains. This is a lattice 
artifact, and it remains a major problem when one studies the restoration of the spon- 
taneously broken chiral symmetry at finite temperature. Moreover, the very successful 
(9(a^) improvement of the pure gauge action makes it very natural to try to bring the 
quark action to the same accuracy. 

This paper is organized as follows: in section 2 we discuss the improvement of both the 
gauge and the fermion actions and the properties of the free fermion actions. In section 
3 we present the results of the simulations and the comparison of the different actions. 
In particular, we study (a) the mN/mp mass ratios at several fixed values of mT^/irtp as 
functions of the lattice spacing, (b) the Lorentz invariance of tt and p meson states, and 
(c) the restoration of the flavor symmetry (as determined by the mass difference of the 
pseudo-Goldstone and non-Goldstone vr mesons) . Our conclusions are presented in section 
4. 

2 Improvement of the action 

2.1 The gauge action 

We generate gauge configurations with the tadpole- improved SU(3) gauge action [^ ||, [12| : 

Sg = Ppl E (1 - ^M-) + /^rt E (1 - ^M-) + /^Pg E (1 - C,ua) (1) 

where P is the standard plaquette in the fi, v -plane, and R and C denote the real part of 
the trace of the ordered product of SU(3) link matrices along 1x2 rectangles and 1x1x1 
paths, respectively: 



pi,v 



R 



■/lU 



ReTr 



-ReTr 



(2) 
(3) 



(-^fiua — 7TJ^6 ^^ 



(4) 



In general, the improvement conditions do not uniquely specify the form of the action. 
For example, at tree-level, adding either the planar 6-link term or one of several 8-link 
terms to the standard action would cancel the O(a^) errors. However, when the quantum 
corrections are calculated with the lattice perturbation theory, then at least two terms 
are required to cancel 0{g'^^a^) errors 0. The terms in Eq. (|I]) provide the most compact 
form of the action. 

Due to the UV divergence of the tadpole-type graphs in lattice perturbation theory, 
operators formally of order a" in the expansion of the action are changed to order a""^"**?^™" 
by quantum effects, depending on the number of tadpole graph contributions to that 
particular term. (In the tadpole contributions n > 2771, i. e. tadpoles do not introduce 
additional UV-divergencies.) The contribution of the tadpole diagrams can be partially 
taken into account by absorbing them in the lattice coupling constants. This is commonly 
achieved by the definition of the 'average gauge link' from the plaquette, Uq = {Py^'^, 
which is strongly dominated by tadpoles, and by replacing f/j(x) -^ Ui{x)/uo in every 
lattice operator 0, [1^. This corresponds to a redefinition of the lattice gauge coupling 
g^ -^ g'^/ul 

With these ingredients, the coefficients of the action (P are related by the 1-loop 
expressions 
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/^rt = --^(1 + 0.4805a,) (5) 

/3pg = _^o.03325a, (6) 

where the strong coupling constant is determined through the 1-loop relation 

a, = -41og(Mo)/3.0684. (7) 

The leading errors of this action are of order 0{a?al, a^). 



2.2 Tree-level improvement of the quark action 

In this work, we study the following fermion action: 

+ C2 p/,(x)?7^(x + /i)t//,(x + 2/i)x(x + 3/i) 

-f/t(x - /i)f/t(x - 2ii)Ul{x - 3^i)x{x - 3/i)] I (8) 

+ a^mgY,xix)xix), 

X 

where the phase factor r]^{x) = (_i)(^o+2:i-- :i>-i)_ xhe standard Kogut-Susskind (stag- 
gered) action is obtained with coefficients Ci = 1 and C2 = 0. At tree level, the action 
is O(a^) accurate when ci = 9/8 and C2 = —1/24. In this case, the difference from the 
Kogut-Susskind action is a discrete version of the 3rd order derivative: 

1 f{x + fi)- fix - ft) 1 f{x + 3/x) - f{x - 3/i) _ a^ ^3 ^^^^ , ^^^^4^ 

8 2a 24 2a ~ 5 '^^^ ^ ^ ^^ ' ^^ 

The staggered action with a third nearest neighbor term was originally proposed by Naik 

[0 . However, he was studying the improvement of the Dirac-Kahler action, which has a 

different coupling to the gauge fields than the action in Eq. (j^). The Dirac-Kahler action 

lacks the exact f/(l) x f/(l) -symmetry enjoyed by the action (|^), and the bare quark mass 

has to be additively renormalized. These properties make use of the Dirac-Kahler action 

much less appealing than the Kogut-Susskind action. Nevertheless, in the following we 

shall call the action (§) the Naik action. 

The (one-component) Grassmann field x describes 4 fiavors of Dirac fermions in the 

continuum limit. This is not transparent in Eq. (^, nor can one easily identify the leading 

irrelevant terms when the continuum limit is taken. At the free fermion level, perhaps the 

easiest way to see this is to use the following transformation [|l^ : in momentum space, we 



decompose the momentum vector k = p+irA/a, A^ = or 1, and we restrict < k < n/a. 
A new (16- component) fermion field ip is defined as 

Hp) = ^E(-l)^''rAX(p + vr5/a) (10) 

^ A,B 

m = ^E(-l)^Mx(p + vri?/a). (11) 



where 

Ta = l^'^lth^'l^' . (12) 

In terms of field ip, the free action (Rl) becomes 



S{ree = J2'^(P) 



^7^- (ci sinp^a + C2 sin 3p^a) + m 



V^(P) • (13) 



This form of the action is flavor diagonal; however, if we perform an inverse Fourier 
transform, the derivative term becomes nonlocal. The Kogut-Susskind action (c2 = 0) 
has leading O(a^) errors. The coefficients ci = 9/8 and C2 = —1/24 for the Naik action 
are readily recovered from Eq. (|13D by expanding the trigonometric functions. 

When the gauge fields are included one cannot transform the action (H) to the form 



in Eq. (13). It is not at all obvious that the interacting Kogut-Susskind action is still 
0{a) -accurate. In order to see the fiavor structure more clearly, one usually performs 
the (local) transformation originally proposed by Kluberg-Stern et al. [Q. It transforms 



the 1-component staggered field x to a hypercubic 16-component 'quark field' (4 fiavors 
of 4-component Dirac spinors), which lives on a lattice with twice the original lattice 
spacing. The quark field action cannot be written in a compact form, but when expanded 
in powers of the lattice spacing a it has apparent dimension-5 terms (giving rise to 0{a) 
errors) . 

However, the Kogut-Susskind action does not have on-shell 0{a) errors. This has been 
shown by Sharpe [|T3| and Luo [|TB], |T7| by a generalization of the Kluberg-Stern et al. 



transformation. The leading scaling violations start at (9(a^). In order to cancel them, 
one has to add dimension-6 terms to the action; these terms have been classified by Luo 
[[17| . The terms fall into two classes: ipD^ip terms, where D^ is a generic combination 
of 3 covariant derivatives (and to which class the 'Naik term' in Eq. (H) belongs), and 
4-fermion terms. Unfortunately, even in the simplest form, the action has 15 dimension-6 
terms with — so far — unknown coefficients. Therefore, we limit ourselves here to a much 
more modest goal and study the degree of improvement possible to obtain with the action 
(P), bearing in mind that this action cannot cancel all of the 0{a'^) errors, but only the 
ones present already for free fermions. 

As with the gauge action, we may improve the action (|^) beyond the tree-level by taking 
into account the modifications due to gluon tadpoles: with the replacement U -^ U/uq, 



the coefficients q in action (|) become 

9 1 

Ci = C2 = ^ . (14) 

In this work we use the quark action defined by Eqs. (|8|.|T^). In the nonzero temperature 
calculation in Ref. |jlO| the action (H) was used with the 'tree-level' values. 



2.3 Properties of the free quark action 

The free quark dispersion relation -E(p) can be found from Eq. ([T3|) , by solving for the 
poles of the Euclidean propagator and using the identification E = Rezpo- In Fig- we 
show the massless quark dispersion relations for the standard Kogut-Susskind and Naik 
actions. For comparison, we also show the Wilson fermion action dispersion relation. 
The Naik action follows the continuum dispersion relation ii^ = |p| much better than the 
standard Kogut-Susskind action up to |p| ~ 1.8/a, not to mention the Wilson action (with 
the Wilson parameter r = 1). Note that for massless free quarks both the Wilson and 
the Kogut-Susskind actions have O(a^) leading errors. Due to the third nearest neighbor 
coupling in the imaginary time direction, unphysical ghost branches (with complex ipo) 
appear in the dispersion relation. These states will become infinitely massive when a — * 0. 

For free fermions, the thermal energy can be easily calculated from E = T'^d log Z/dT. 
In the imaginary time formalism, the temperature T = l/{NTa), the inverse temporal 
extent of the lattice. Under the assumption that the free energy is proportional to the 
volume, the pressure is P = Td log Z/dV. In Fig. ^ we show E/T^ and P/T^ for free 
Kogut-Susskind, Wilson and Naik fermions as functions of the inverse lattice spacing. 



Also shown are the results from the Bielefeld "P4" staggered action ||T8|: like the Naik 
action, it contains a 3rd nearest neighbor coupling, but in this case the neighbors are 
coupled along L-shaped paths (the Naik 'ijj{x)U^{x)U^{x + fi)U^{x + 2}l)ip{x + 3/1) -terms 
are replaced with terms of form iIj{x)U^{x)U^{x + fi)U^{x + fi + i))ip{x + (1 + 2z>), with 
u y^ fi). The P4 action yields the same tree- level improvement as the Naik action. (For a 
comparison with a renormalization group improved free staggered quark action, see the 
last paper of Ref. 0.) 

The energy and pressure of the Naik fermions approach the continuum ideal fermion gas 
limits much faster than the standard Kogut-Susskind action. Indeed, the Bielefeld group 
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Figure 1: The dispersion relation £^(p) for massless free quarks with different fermion 
actions. The momentum p is to the spatial direction (1, 1, 0), and the dispersion relations 
are plotted up to the end of the Brillouin zone. 



[|T0] reported an improved thermodynamic behavior even when the interacting gauge fields 
are included in a dynamical quark Monte Carlo simulation. 



3 The simulations and the results 

3.1 Hadron spectrum 

The parameters of the action (|l]) used in the generation of the quenched configurations 
are shown in Table |l|. 



/3pi Up 



volume A*", 



conf. 



6.8 0.8261 16^ X 32 199 
7.1 0.8441 143 X 28 203 
7.4 0.8629 16^ x 32 200 
7.6 0.8736 16^ x 32 100 
7.75 0.8800 16^ x 32 200 

7.9 0.8848 16^ x 32 200 



Table 1: The parameters of the runs. (5-ct and /5pg can be obtained through Eqs. (§-0) 



t3 



0.3 




\ Wilson 



Kogut-Susskind 



Naik 
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12 16 
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0.0 I • ^ 

24 4 



12 16 
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Figure 2: The energy (left) and the pressure (right) per fermion degree of freedom for 
free Kogut-Susskind, Naik, Wilson and "P4" |T^ fermions as a function of Nt = l/{aT). 
The continuum values are shown with dashed lines. 



We measure the masses of the nucleon, the Goldstone pion tt, (corresponding to the 
(spontaneously broken) explicit U(l) chiral symmetry of the action (|^)), the non-Goldstone 
("SC") pion 7i2, and the p and p2 mesons. The masses were calculated both with the Naik 
and the Kogut-Susskind actions. For each lattice and propagator, we use four wall source 
planes. In each case, the hadron propagators were measured with 5-6 bare quark masses 
aruq = 0.005 - 0.32; the hadron masses are shown in Tables |^-^. 

Throughout the analysis we quantify the performance of the improved actions by com- 
paring the results against a non-improved benchmark — an extensive standard quenched 
Kogut-Susskind hadron spectroscopy study by the MILC collaboration [|lT|. In particular, 
we use /?wiison = 6/^^ = 5 54 (^gS)^ 57 (24^), 5.85 (24^) and 6.15 (32^) lattices (with the 
spatial volume in parentheses). 

The Naik hadron propagator calculation requires about 2 times more CPU time than 
the Kogut-Susskind one. The number of conjugate gradient iterations is very similar for 
the Naik and the Kogut-Susskind quarks, but since the Naik action (|^) involves about 
twice as many terms, the computational load is higher. For example, for the /5pi = 7.4, 
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Naik, 


/3pi = 6.8, 16^ X 32 








aiTLq 


71 1X2 


P 


P2 


Nucleoli 


0.02 


0.343438(70) 1.4(2) 


1.537(15) 


2.06(18) 


2.301(50) 


0.04 


0.484458(64) 1.53(3) 


1.617(24) 


2.021(72) 


2.387(19) 


0.08 


0.682699(82) 1.86(6) 


1.6964(91) 


2.098(29) 


2.730(38) 


0.16 


0.962771(60) 2.12(3) 


1.8581(75) 


- 


3.032(49) 


0.32 


1.365621(51) - 


2.061(27) 


- 


3.589(14) 



Kogut-Susskind, /3pi = 6.8, 16^ x 32 



arriq 


vr 7l2 


P 


P2 Nucleoli 


0.02 


0.350547(81) 1.411(75) 


1.405(13) 


1.73(13) 2.1632(93) 


0.04 


0.492871(80) 1.474(35) 


1.4308(61) 


1.642(48) 2.134(64) 


0.08 


0.689935(67) 1.524(16) 


1.4816(25) 


1.660(17) 2.3286(59) 


0.16 


0.959548(63) 1.640(10) 


1.5772(37) 


1.7389(56) 2.447(33) 


0.32 


1.325654(74) 1.8896(43) 


1.75802(52) 


1.9078(25) 2.7434(79) 



Naik, 


f3pi = 7.1, 143 X 28 








arrig 


vr 7r2 


P 


P2 


Nucleon 


0.02 


0.35637(12) 1.289(39) 


1.573(44) 


1.595(67) 


2.16(13) 


0.04 


0.50074(20) 1.473(59) 


1.568(18) 


1.651(38) 


2.342(51) 


0.08 


0.70269(17) 1.585(26) 


1.6188(70) 


1.822(23) 


2.459(31) 


0.16 


0.985219(93) 1.968(65) 


1.822(12) 


2.0087(83) 


2.759(19) 


0.32 


1.389359(81) 2.874(80) 


2.1188(39) 


2.72(16) 


3.384(74) 



Kogut-Susskind, /3pi = 7.1, 14'^ x 28 



arug 


TT 


7r2 


P 


P2 Nucleon 


0.02 


0.36368(13) 


1.411(77) 


1.338(14) 


1.376(31) 1.990(37) 


0.04 


0.50958(13) 


1.356(32) 


1.403(18) 


1.443(16) 2.190(59) 


0.08 


0.70934(13) 


1.361(25) 


1.4481(69) 


1.5301(72) 2.278(18) 


0.16 


0.97939(11) 


1.592(17) 


1.5680(76) 


1.6681(53) 2.414(11) 


0.32 


1.34259(10) 


1.8213(89) 


1.7423(15) 


1.8574(16) 2.7142(25) 



Table 2: Masses for /3pi = 6.8, 16^ x 32, and f3pi 
no good mass fits were possible. 



7.1, 14^ X 28 lattices. Entry '-' means 



16"^ X 32 lattices the number of the conjugate gradient iterations for each source plane 
varies approximately from 130 {arrig = 0.32) to 2050 {arrig = 0.02) for the Kogut-Susskind 
and from 140 to 2400 for the Naik action, whereas the CPU time per plane for arug = 0.02 
is about 260 seconds for K-S and 600 seconds for Naik on the Intel Paragon using 32 nodes. 
In order to find the best confidence levels of the propagator fits, we used one-, two- and 
three-particle fitting functions, varying both the beginning and the end of the fitting range. 
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Naik, 


/3pi = 7.4, 163 


X 32 




aiTLq 


TT 


71-2 p 


P2 Nucleoli 


0.02 


0.37314(79) 


1.033(30) 1.268(13) 


1.294(34) 1.765(30) 


0.04 


0.52238(70) 


1.0834(96) 1.3139(89) 


1.383(14) 1.920(15) 


0.08 


0.72475(95) 


1.2517(75) 1.4283(99) 


1.5035(71) 2.1320(68) 


0.16 


1.0080(11) 


1.5591(97) 1.6241(48) 


1.768(14) 2.547(14) 


0.32 


1.416801(79) 


2.146(23) 1.9719(25) 


2.2396(79) 3.1459(76) 



Kogut-Susskind, /3pi = 7.4, 16^ x 32 



arriq 


vr 


^2 P 


P2 Nucleoli 


0.02 


0.38080(74) 


0.9603(66) 1.207(14) 


1.247(33) 1.724(29) 


0.04 


0.5280(12) 


1.0418(45) 1.2449(63) 


1.337(43) 1.855(13) 


0.08 


0.7297(10) 


1.1905(66) 1.3270(28) 


1.4074(63) 2.0282(53) 


0.16 


0.9994(11) 


1.4197(28) 1.4774(25) 


1.5725(50) 2.3174(56) 


0.32 


1.362520(74) 


1.7300(16) 1.69637(81) 


1.7975(31) 2.6656(35) 



Naik, /3pi = 7.6, 16^ x 32 



ania 



0.01 
0.02 
0.04 
0.08 
0.16 
0.32 



vr 



vr2 



P 



P2 



Nucleoli 



0.27294(23) - 0.945(28) - 1.36(15) 

0.38154(24) 0.7513(89) 1.0106(87) 1.046(11) 1.411(52) 

0.53125(25) 0.8629(46) 1.0911(56) 1.1424(73) 1.598(11) 

0.73613(23) 1.0559(32) 1.2420(53) 1.312(11) 1.852(40) 

1.01922(24) 1.3667(35) 1.4682(31) 1.5441(56) 2.256(15) 

1.42502(16) 1.9038(93) 1.8670(34) 2.0072(69) 2.9137(61) 



Kogut-Susskind, /3pi = 7.6, 16^ x 32 



arrig 


TT 


71-2 p 


P2 Nucleon 


0.01 


0.28090(23) 


0.945(29) 


1.28(14) 


0.02 


0.39165(22) 


0.7508(90) 0.9943(86) 


1.031(11) 1.377(54) 


0.04 


0.54237(23) 


0.8642(48) 1.0621(48) 


1.190(15) 1.577(10) 


0.08 


0.74470(22) 


1.0462(31) 1.2030(44) 


1.2705(91) 1.804(33) 


0.16 


1.01321(19) 


1.3092(27) 1.3882(22) 


1.4474(38) 2.187(20) 


0.32 


1.36891(12) 


1.6677(32) 1.6573(10) 


1.7287(19) 2.5975(52) 



Table 3: Masses for /3pi = 7.4, 14^ x 28, and /3pi = 7.6, 16^ x 32 lattices. 

All of the fits use the full invariance matrix of the propagators. We block together all of 
the propagators on each lattice, then, in order to facilitate further analysis, we calculate 
the masses using a single elimination jackknife procedure. When fitting each jackknife 
sample, we use the invariance matrix of the entire ensemble, rather than recomputing the 
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Naik, /3pi = 7.75, 16^ x 32 



aniq 


TT T^2 P P2 


Nucleoli 


0.01 


0.26943(34) 0.5374(51) 0.829(18) 0.862(12) 


1.140(80) 


0.02 


0.37638(30) 0.6175(42) 0.8951(96) 0.906(11) 


1.265(26) 


0.04 


0.52339(20) 0.7421(22) 0.9701(47) 1.015(13) 


1.429(10) 


0.08 


0.72718(18) 0.9428(14) 1.1071(22) 1.1479(77) 


1.6855(35) 


0.16 


1.01186(16) 1.2587(17) 1.3577(65) 1.3963(26) 


2.1018(42) 



Kogut-Susskind, /jpi = 7.75, 16^^ x 32 



arriq 


vr 7^2 P P2 Nucleoli 


0.01 


0.27794(33) 0.5445(50) 0.830(16) 0.863(11) 1.174(11) 


0.02 


0.38752(28) 0.6182(22) 0.8895(82) 0.8898(53) 1.282(10) 


0.04 


0.53656(21) 0.7504(21) 0.9632(41) 1.025(11) 1.4225(90) 


0.08 


0.73830(17) 0.9460(13) 1.0984(32) 1.1351(54) 1.6866(64) 


0.16 


1.00860(12) 1.2320(16) 1.3173(32) 1.3550(41) 2.0571(47) 



Naik, /5pi = 7.9, 16^ x 32 



arUn 



0.005 

0.01 

0.02 

0.04 

0.08 

0.16 



71 



vr2 



P 



P2 



Nucleon 



0.18486(32) 
0.25907(28) 
0.36229(28) 
0.50589(26) 
0.70786(24) 
0.99501(20) 



0.3806(67) 
0.4284(48) 
0.5054(37) 
0.6400(20) 
0.8506(16) 
1.1710(14) 



0.6681(99) 

0.760(23) 

0.7873(98) 

0.8655(89) 

1.0056(33) 

1.2646(32) 



0.6831(75) 

0.725(11) 

0.7678(53) 

0.8546(49) 

1.0104(28) 

1.2790(17) 



0.920(12) 

1.017(16) 

1.0972(75) 

1.2632(72) 

1.5236(93) 

1.9560(52) 



Kogut-Susskind, /3pi = 7.9, 16^ x 32 



aiTiq 


TT T^2 P P2 Nucleon 


0.005 


0.19160(32) 0.3848(71) 0.6704(96) 0.6756(73) 0.9114(67) 


0.01 


0.26809(28) 0.4339(49) 0.720(24) 0.753(28) 1.005(15) 


0.02 


0.37383(31) 0.5142(36) 0.7767(87) 0.796(11) 1.0963(69) 


0.04 


0.51985(29) 0.6525(20) 0.8599(44) 0.8576(45) 1.2703(65) 


0.08 


0.72125(26) 0.8624(15) 1.0051(23) 1.0115(26) 1.5376(78) 


0.16 


0.99577(21) 1.1642(27) 1.2401(19) 1.2630(37) 1.9388(54) 



Table 4: Masses for /3pi = 7.75 and l3^\ = 7.9, 16^ x 32 lattices. 

invariance matrix for each sample. 

A generic feature of the fits to the propagators is that one has to use considerably larger 
minimum fit distance from the source with Naik fermions than with the Kogut-Susskind 
fermions. The reason for this effect is probably the large extent in the imaginary time 
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Figure 3: Nucleon (upper) and p (lower) masses as functions of aruq for (3^\ = 7.4 and 7.9. 

16^x32 )Spi=7.4 
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Figure 4: The pion mass squared as a function of aniq for /3pi = 7.4. 

direction of the Naik derivative operator in Eq. (§) [|l^]. The transfer matrix is well 
defined only at imaginary time separations larger or equal to 3. The ghost branch in the 
dispersion relation can also cause short distance effects in the correlation function. 

In Fig. 1^ we summarise the nucleon and p meson masses from Tables ^ and | for 
/3pi = 7.4 and 7.9, and in Fig. | the pion mass squared for /?pi = 7.4. The masses of the 
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Figure 5: The Edinburgh plots for (3p\ = 7.4, 7.6, 7.75 and 7.9. The fancy crosses show 
the standard unimproved Kogut-Susskind data at /3wiison = Q/g^ = 5.54, 5.7 and 5.85 | PT| ; 
these correspond roughly to the same lattice spacing (determined by arrip) as the improved 
/3pi = 7.4, 7.75 and 7.9. The small circles denote the physical limit {niT^/mp ~ 0.18) and 
the infinite quark mass limit {m-,^/mp = 1). 

Naik hadrons in lattice units tend to be larger than the Kogut-Susskind masses, but the 
difference gets smaller with decreasing arriq (approaching the chiral limit) and increasing 
/3pi (continuum limit). At /5pi = 7.9 the differences are barely discernible. 
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Figure |^ shows the Edinburgh plots for /3pi = 7.4-7.9. The m^ jnip -ratios from /3pi = 
6.8 and 7.1 exhibit typical strong coupling behavior: the ratio vtij^jvtip remains roughly 
constant at around 1.5 when vrtq -^ [711^^/ nip -^ 0). Only when (3p\ > 7.4 do the 
data show an approach to the vicinity of the physical value, and thus in the following 
we concentrate on these couplings. For comparison, we also plot the m-jv/w^p mass ratio 
obtained with standard non-improved Kogut-Susskind action at P = 5.54, 5.7 and 5.85 
[|lT) ; these couplings correspond roughly to the same lattice spacing as the improved action 
at /3pi = 7.4, 7.75 and 7.9. 

It is interesting to note that despite the large differences in the masses in lattice units 
in Fig. I (and Table 0), in Fig. |^ the /?pi = 7.4 Naik and Kogut-Susskind mass ratios lie 
practically on the same curve, with the Naik values displaced slightly in the direction of 
smaller niT^/mp. When the mass of the p meson is used to set the scale, from Figs. ^ and 
^ we see that for a given set of bare parameters the lattice spacing for the Naik fermions 
is slightly larger than for the Kogut-Susskind fermions, while the mass ratios are closer to 
the physical values. Conversely, if we want to investigate the same physical system (mass 
ratio and physical volume), the Naik action enables us to use a bit larger bare quark mass 
and smaller lattices (in lattice units). This effect becomes smaller as one gets closer to 
the chiral limit irtq —>■ and to the continuum limit /3pi -^ 00; nevertheless, it is clearly 
observable in all of the Edinburgh plots in Fig. ^. 

3.2 The chiral function and the continuum Umit 

A convenient method to quantitatively measure the degree of improvement in hadron 
spectroscopy is to study the lattice spacing dependence (in units of anip) of the ratio 
mi\i/mp at some fixed value of m.,^/mp. This requires interpolation or extrapolation to the 
desired m-^/mp ratio. We perform this for each f3p\ separately with chiral fit functions. 

The chiral fits are motivated by quenched chiral perturbation theory {QxPT), which 
gives rriN and rup in a power series of m.,^ (+ logarithmic terms). Since in the leading 
order m^ oc nig, these become power series in m^"^. For the standard Kogut-Susskind 
action, the chiral extrapolations have been discussed in detail in [|ll|]. 

Extrapolation of the ratio mi^/mp to the chiral limit m^ — i> or even to the physical 
limit iriT^/nip ^ 0.1753 is sensitive to the form of the selected chiral fit function ansatz. 
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The value of the ratio is much less sensitive in the region mT^/nip ^ 0.4-0.6, where the 
function interpolates between measured mass values (see Fig. ^). Detailed comparison of 
the different actions is feasible in this region. 

We exclude the "strong coupling" runs at /3pi = 6.8 and 7.1, and fit amN^amq) and 
amp{amq) for /Spi = 7.4-7.9 to the chiral ansatz 

am = Co + Ci aruq + C3/2(amg)^'^ + C2{amqY . (15) 

This function gives good fits at all 4 couplings (after excluding the anomalous smallest 
arriq value 0.005 from /3pi = 7.9). 

QxPT for nucleons and vector mesons ||2^ implies the presence of an additional term 
oc am^, which corresponds to {amqY^'^ in the continuum. If one includes a term oc (arrig)^/^ 
in ([XSp, the fits invariably prefer a positive sign for the coefficient; whereas QxPT gives a 
negative sign. However, the appropriate pion mass in this term is actually the fiavor singlet 
pion mass, which is not proportional to (ar/ig)^/^ at fixed lattice spacing due to fiavor 
symmetry breaking. When this is taken into account, acceptable fits with a coefficient 
compatible both in sign and magnitude with QxPT are possible. This is studied in detail 
for our standard gauge Kogut-Susskind hadrons in Ref. [|lT|. Since such fits do not appear 



to change the extrapolated/interpolated values significantly from (|T^), but do increase the 
errors, we prefer to leave out the amj^ term. 

The error propagation is taken into account by performing the (fully correlated) fits 
separately to each of the jackknife blocks. 

In most cases, it would be possible to obtain acceptable fits also with the simpler ansatz 



with either C3/2 or C2 fixed to zero. However, while the full ansatz ([TsD works quite well 



for the standard Kogut-Susskind hadrons [|r^, these simplified functions do not. In order 



to facilitate the comparisons between the different actions, we retain the full chiral ansatz 
(|15D here. 

The results of the chiral extrapolation/interpolation to 771^^/ nip = 0.1753 (physical), 
0.4, 0.5 and 0.7 are shown in Fig. P; both for improved actions and for the standard 
Kogut-Susskind action. Since we expect the leading errors to be C(a^), we plot the ratios 
against {ampY'. Here airip is calculated at the quark mass which yields the indicated value 
of mT^/irip. We make a linear fit with respect to {anipY of the standard Kogut-Susskind 
data, and, since in the continuum limit all of the actions must yield equivalent results, 
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Figure 6: The itln /nip mass ratios as functions of the square of the lattice spacing (in 
units of {anipY), for rriT^/mp = 0.1753, 0.4, 0.5 and 0.7. From left to right, the Naik and 
the (improved gauge) Kogut-Susskind points correspond to Pp\ = 7.9, 7.75, 7.6 and 7.4; 
the standard Kogut-Susskind points to /Swuson = 6.15, 5.85, 5.7 and 5.54. The straight 
lines are linear fits to the (from top to bottom) standard K-S, improved K-S and Naik 
data, where the intercept at a = in the last two fits is fixed to the standard K-S value. 

we fit straight lines to the improved Kogut-Susskind and Naik data, with the constraint 
that the a = intercept is fixed to the standard Kogut-Susskind value. 
We make the following observations: 
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In the intermediate m.^/mp = 0.4 and 0.5 plots, the improved gauge nucleon to p 
mass ratios are clearly closer to the continuum values than the standard Kogut- 
Susskind ones. Indeed, the /3pi = 7.9 value is very close to the standard Kogut- 
Susskind /Swuson = 6.15 one, but with twice the lattice spacing (albeit with larger 
statistical errors). At large lattice spacings {jSpi = 7.4) the Naik fermions show 
smaller scaling violation than the improved gauge Kogut-Susskind fermions, but 
this difference becomes very small when the lattice spacing is reduced. 

At the physical ratio m-j^/mp = 0.1753 the errors in nii^/mp increase dramatically 
due to the extrapolation in amq. Nevertheless, we observe a pattern similar to that 
at larger quark mass. 

When iriT^/mp = 0.7 the quark mass arug becomes so large that the chiral expansion 
( [T3| ) does not converge well any more: the highest power terms have the largest 
magnitude. This leads to erratic jumping of the points in the last panel of Fig. p. 
(note however the very small range of m^ /irip covered by this plot). 

As the quark mass is lowered, the difference between the two types of quarks in the 
improved gluonic fields is reduced. Thus, at m.,^/mp = 0.7 most of the improvement 
of niN/iTip is attributable to the Naik improvement, whereas near the physical quark 
mass, most of the improvement comes from the gluonic action. A large part of the 
Naik improvement is due to the larger (amp) and hence a larger lattice spacing. If 
one uses the string tension to set the scale the difference between the Naik and the 
Kogut-Susskind actions becomes smaller. 

The linearity (against {ampY) of the standard Kogut-Susskind mi^/mp-idXio clearly 
supports the notion that the scaling violations behave as 0{a?). When mT^/rrip < 
0.5, the constrained linear fits to the improved gauge Kogut-Susskind and Naik 
data have confidence levels better than 0.5, certainly quite compatible with O(a^) 
leading scaling violations. The magnitude of the violations - the slope of the line - 
for the Naik data is only about 1/2 of the standard Kogut-Susskind value, whereas 
the improved gauge Kogut-Susskind has a slightly larger slope than Naik. 

We can also test whether the data would allow for O(a^) scahng violations for the 
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Naik action. When mT^/nip = 0.5 a constrained fit of form A + B{amp)^ (where 
again A is set to the a = intercept of the standard Kogut-Susskind data) does not 
fit the Naik data well: the confidence level is only 0.15 (as opposed to 0.75 before). 
This disfavors the leading 0{a^) errors. For smaller 771^^/ rrip-T atios the statistical 
errors become larger and this analysis cannot distinguish the fits. 

• To check consistency, we can relax the constraint at a = and fit independent 
straight lines to all datasets. When mT^/nip < 0.5 the intercepts at a = are 
compatible for all cases, i. e. , within 1 standard deviation of each other. 

The lattice spacing and the size of the system in physical units can be obtained by 
extrapolating anip to the physical mT^/irip-iaXio and setting mp = 770 MeV. These are 
given in Table ^. 



Ppl 


6.8 


7.1 


7.4 


7.6 


7.75 


7.9 


a (fm) 
Size (fm) 


0.37 
5.9 


0.35 
5.0 


0.31 
5.0 


0.24 
3.8 


0.21 
3.4 


0.19 
3.0 



Table 5: The lattice spacing and the box size in physical units. 

The numbers in Table ^ have been calculated with the Naik quark action; for the Kogut- 
Susskind action the lattice spacings and the box sizes would be fractionally smaller. The 
box sizes are considerably larger than 2fm (with the possible exception of (3p\ = 7.9), so 
that we can safely ignore the finite size effects ||2I|. For the weakest coupling and the 
smallest quark mass, the product m^ x (Lattice size) is approximately 3.0. 

Besides the mass of the p-meson, the square root of the string tension is commonly 
used to determine the lattice spacing. In Fig. |^ the ratio nip/y/a is plotted against a^a. 
Here arrip is evaluated at the physical mT^/nip = 0.1753 and at 0.7. The string tension 
for the standard gauge action is interpolated from the data in the literature |2^; for the 
improved gauge action (|l|) it has been measured by the SCRI group |2^ . 



Since the scale violations are expected to behave as 0{a ), the rrip/y/a ratio should 
behave linearly as a function of a'^a. Indeed, the standard Kogut-Susskind data shows 
good linearity up to a^a = 0.17 (/3wiison = 5.7). However, at stronger coupling (5.54) 
the ratio falls strongly off the linear behavior. We fit a straight line to the three weakest 
coupling datapoints, the intercepts at a = are 1.738(25) at the physical m-^/mp = 
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Figure 7: The ratio nip/ y/a against a?a at mT,/mp = 0.1753 and 0.7. 



0.1753 and 2.238(11) at 0.7 (the errors quoted here are only statistical). These results are 
consistent with the SCRI group (preliminary) Wilson and clover fermion mass ratios p3| . 

The improved gauge Kogut-Susskind and Naik data still seem to reside completely in 
the "strong coupling region" , although there is some indication that at weaker couplings 
the ratios would bend to the direction of the line defined by the standard Kogut-Susskind 
data. Extrapolation of the improved action ratios to the continuum limit is clearly not 
justified. 

The non-linearity in vtipj ^ is somewhat surprising, when we compare it against the 
purely hadronic observables in Fig. ^. This lends support to the view that a large part 
of the scaling violations cancel in the hadronic ratios, and justifies the use of arrip as the 
scale factor in purely hadronic observables. 



3.3 Lorentz symmetry 



As discussed in Sec. ^]^, the free quark continuum dispersion relation is approximated 
much better by the Naik action than by the standard Kogut-Susskind action. At very high 
temperatures, deep in the quark-gluon plasma phase, the quarks behave approximately 
as free particles, and the Naik action improves QCD thermodynamics [l^. However, a 
priori it is not clear whether the dispersion relation of hadronic states is improved. 
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Kogut-Susskind vr 


Naik 7r 


kL/(27r) 


aE{k) c2(k) 


aE{k) 


c\k) 


(0,0,0) 


0.53521(17) - 


0.52625(15) 


- 


(0,0,1) 


0.65223(56) 0.9010(45) 


0.65111(50) 


0.9532(39) 


(0,1,1) 


0.74710(63) 0.8809(30) 


0.75474(76) 


0.9489(36) 


(1,1,1) 


0.82655(88) 0.8575(31) 


0.8411(11) 


0.9305(40) 


(0,0,2) 


0.8817(14) 0.7959(40) 


0.9087(26) 


0.8897(77) 


(0,2,2) 


1.0925(32) 0.7353(57) 


1.1549(44) 


0.8566(83) 


(2,2,2) 


1.275(11) 0.723(15) 


1.395(19) 


0.902(29) 





Kogut-Susskind p 


Naik p 


kL/(27r) 


aE{k) 


c\k) 


aE(k) 


c\k) 


(0,0,0) 


1.2411(69) 


- 


1.3065(79) 


- 


(0,0,1) 


1.262(27) 


0.34(44) 


1.3489(80) 


0.73(15) 


(0,1,0) 


1.289(13) 


0.79(21) 


1.362(16) 


0.97(28) 


(0,1,1) 


1.298(13) 


0.47(11) 


1.385(19) 


0.68(16) 


(1,1,0) 


1.320(18) 


0.65(16) 


1.385(12) 


0.68(12) 



Table 6: The energy of the tt and p meson states at finite momentum k = wl-n jL^ and the 
'speed of hght squared' c^i^) = {E^k) - E\0))/k'^, for /3pi = 7.4, 



arrin 



0.04, 16^^ X 32 



lattice. 



We test hadron dispersion relations by measuring the energy of the tt- and p-meson 
states with finite spatial momenta on 100 lattices with f3p\ = 7.4, aruq = 0.04 and volume 
16^ X 32. We use 4 (finite momentum) wall sources, separated by 8 lattice units. 

The source operators are constructed as follows: first, we take a zero momentum wall 
source, which is 1 for a particular source color at each spatial slice at the source time. 
This is used as a source for the conjugate gradient to compute the quark propagators. 
Then this wall source is multiplied by the momentum dependent phase factor exp(ik • x), 
by the sign factors (depending on the location in the 2*^ flavor hypercube) to select the 
desired meson, and by an extra (— l)^**^^ corresponding to 75. This is used as a source 
for the conjugate gradient to compute the antiquark propagator. The sink operator is 
similar, except that the quark and antiquark propagators are multiplied together with 
the appropriate phase and sign factors before summing over spatial points, corresponding 
to a local sink. 

For pions, we use momentum vectors pointing to 3 different directions: kL/(27r) = 
(0,0, 1), (0,1,1), (1,1,1), and these multiplied by 2. For the p-meson, we use (the lattice 
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Figure 8: The 'speed of light squared', calculated from the pion dispersion relation, for 
Naik and K-S pions. 

analog of) the vector operator ip^sip, and we expect that the dispersion relation may 
be different along the z-axis direction and perpendicular to it. Therefore, for p we use 
kL/(27r) = (0,0, 1), (0,1,0), (0,1,1) and (1,1,0). The signals for higher momenta are too 
noisy to be useful. The results are listed in Table ^. 

The violation of Lorentz invariance can be quantified by measuring the 'speed of light' 
with the continuum dispersion relation 



c\k) 



k2 



(16) 



The deviation of c^ from unity directly measures the violation of Lorentz invariance. The 
results are shown in Table ^ and in Fig. ^ (for pions). The Naik pions show a clear 
improvement of c^ over the Kogut-Susskind pions: the deviation from unity is reduced 
approximately by half. The results for the p-mesons seem to indicate a dependence on 
the direction of the momentum (parallel or perpendicular to z). Also, here the c^ is closer 
to unity for the Naik mesons; however, the statistical errors are so large that we cannot 
make definite statements about the improvement. 
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Figure 9: Flavor symmetry breaking parameter 5^ = (m^^ — rn^) / {m'l — m'^) , interpolated 
to iriT^/mp = 0.5, as a function of {amp)"^. The data correspond to the same values of f3 as 
in Fig. ^. The straight lines are linear fits to 2 (standard gauge) or 3 (improved gauge) 
points with the smallest lattice spacings, constrained to go through the origin. 



3.4 Flavor symmetry 

The restoration of flavor symmetry can be discerned by investigating the mass differences 
between vr and tt2 mesons. The first particle is the Goldstone boson corresponding to the 
spontaneously broken U(l) xU(l) chiral symmetry and it becomes massless when arriq — ^ 
even at a finite lattice spacing (Fig. ^). In comparison, the tt2 mesons remain massive in 
the chiral limit, and become massless only when both the chiral and the continuum limits 
are taken. 

We use the dimensionless quantity 



2 2 

nii — nii 



(17) 



to measure flavor symmetry breaking. For the standard Kogut-Susskind quark action, 
this quantity is almost independent of the bare quark mass airtq at small lattice spacings. 
In Fig. 1^, we show 5^ for /5pi > 7.4 improved gauge Naik and Kogut-Susskind hadrons. 
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together with the unimproved glue Kogut-Susskind values, as functions of {ampY . The 
data is interpolated to mT,/mp = 0.5 (compare to the third panel in Fig. |). With this 
constraint, the flavor symmetry breaking parameter reduces to 5.,^ = (^^^2/"^^ ~ l)/3- 

At this value of uiT^/mp, we observe that the flavor symmetry violation at small a is 
reduced by ^ 45% due to the improved gauge action. When the Naik fermions are used, 
(^TT is slightly smaller than with the Kogut-Susskind fermions. However, this situation 
would become reversed, if we used ay/a instead of amp to set the scale. 

Figure ^ clearly indicates that the leading flavor symmetry breaking terms are propor- 
tional to a^ for all of the actions studied. The region linear in {anipY extends to larger 
lattice spacings with the improved gauge. 

A successful additional improvement of the Kogut-Susskind flavor symmetry is the 
MILC "fat link" fermion action [0. That action substitutes the standard gauge links 
with smeared average links in the fermion hopping terms. The averaging process im- 
proves the flavor symmetry dramatically, the improvement being roughly comparable 
both for the standard Kogut-Susskind and the Naik action. These observations indicate 
the importance of the coupling of fermions to the gauge fields for the flavor symmetry of 
the staggered action. The Naik action (j^) can be interpreted naively as a straightforward 
improvement of the (free) fermion dispersion relation. 

The improvement in flavor symmetry from the Symanzik improved gauge action, like 
the improvement from the fat link quark action, can be understood as a suppression of the 
effects of high momentum gluons. Gluons with momentum near vr/a scatter quarks from 
one corner of the Brillouin zone to another, which is roughly equivalent to changing their 



flavor [^, |2^ . The suppression of the high momentum gluons becomes evident when the 
gauge action is expanded to quadratic order in the vector potential A^ (where the lattice 
variable U^{x) = exp[—igaAp_{x)] ). Using the shorthand notations 

A;^ = 2sin-aA;^, P = ^ A;^ , (18) 

and 

UAk) = kMk) - KApik) , (19) 



the quadratic part of the action (|l|) can be written in the form p6| 
^(2) = la2 ^ f^Ak)fpA-k) [cpi + 8crt + 16cpg - (Crt - Cpg)(A:J + kl) - CpgPj . (20) 



k;fi<v 
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Here the coefficients q denote the relative strength of the three terms in the action. They 
are related to coefficients /?j through Ci'ojg^ = (3i. As an overall normalization we require 
that the constant term within the brackets equals to one: Cpi + Scrt + 16cpg = 1. 

For the Wilson gauge action the coefficients are Cpi = 1, Crt = Cpg = 0, whereas for the 
improved action Cpi > 1 and c^t, Cpg < 0. With the improved gauge, the non-constant 
terms within the brackets in Eq. ( pOD increase the action for modes close to the edges of 
the Brillouin zone (fc^ ~ ±71 /a, for at least one /i). 

As a simple example we consider a momentum vector parallel to one of the lattice axes 
and at the edge of the zone. In this case the term within the brackets in Eq. ([20| ) reduces 
to [1 — 4crt]. At /5pi = 7.4, using Eqs. (|,§|), Table |1| and the normalization condition 
above, the coefficent Crt has a value ~ —0.26. When compared to the Wilson gauge 
action (crt = 0), this more than doubles the action difference of the modes close to the 
edge of the zone and near the origin k = 0. When the lattice spacing is reduced the 
suppression of the modes near the edge of the zone increases rapidly, while the relative 
difference between the actions becomes smaller. At tree level, the coefficients assume 
values Cpi = 5/3, c^t = —1/12 and Cpg = 0, still yielding a 33% difference of the action at 
k^ = n/a. 

4 Conclusions 

We investigate improvement of the quenched light hadron mass spectrum using a tadpole- 
improved staggered Naik action (P), which at the tree level does not have O(a^) errors. 
Correspondingly, we use C(a^) tadpole improved gauge action to generate the gauge con- 
figurations. Using the same gauge action for both the Kogut-Susskind and the Naik 
calculations allows us to separate the effect of the fermionic improvement from the im- 
provement of the gauge action. The latter is studied by comparing the results presented 
here to our standard Kogut-Susskind results [[Til] . 

We find that improvement of the gauge action has a significant effect on the hadron 
spectrum: when m.^/mp ~ 0.5, the nucleon to p-meson mass ratio is roughly 50 % closer to 
the continuum value with the improved gauge than with the standard gauge action. Thus, 
the scale violations with the standard gauge spectroscopy are at roughly the same level 
as with the improved gauge at about 1.4 times the lattice spacing. Using the improved 
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gauge action, the Naik quark action has smaller scaling violations than the Kogut-Susskind 
action, although the difference becomes small when the quark mass is reduced. Similarly, 
improving the gauge action reduces the amount of flavor symmetry breaking, but using 
the Naik action yields little further gains. For both of the actions the flavor symmetry 
can be further improved with the 'fat link' procedure 0. 

The biggest improvement provided by the Naik action comes from the improved Lorentz 
invariance of the hadronic states. This is best evidenced by the vr-meson dispersion 
relation, which is much closer to the continuum behavior when the Naik action is used. 
This property may be especially significant for nonzero temperature simulations, where 
the hadronic and/or quark degrees of freedom typically have large momenta. Thus, when 
one strives for higher precision in staggered quark simulations, an economical solution can 
be found from the combination of an improved Naik-like quark action together with the 
fat links. 
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